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1 Introduction 



The existence of a new type of matter, non-baryonic dark matter, is clearly indicated by 
observations on small and large astronomical scales. In particular, the observed flat ro- 
tation curves in spiral galaxies, CMB anisotropics, and matter power spectrum suggest 
the existence of a cold (or warm) dark matter component with energy density around 
five times larger than ordinary baryons[l]. These inferences aside, the precise nature 
of dark matter is not yet known. If dark matter interacts with ordinary matter then it 
might be directly detectable with Earth based experiments, some of which have yielded 
some impressive results. Most notable, are the DAMA/Nal [2] and DAMA/LIBRA [3] 
experiments. These experiments have observed an annual modulation in the 'single hit' 
event rate, over more than 12 annual cycles, with phase consistent with expectations 
from dark matter interactions [4]. Over the last few years, the CoGeNT experiment [5, 6], 
using a p-type point-contact Germanium detector with very low energy threshold, has 
observed a rising event rate at low energies. These events could not be explained 
by known backgrounds and provide evidence supporting the DAMA signal. Most re- 
cently, the CRESST-II collaboration[7] have obtained results that cannot be explained 
by known backgrounds and are compatible with a dark matter signal also rising at low 
energies. 

A satisfactory explanation of these experiments requires a specific model for dark 
matter. One promising approach appears to be in the framework of multi-component 
hidden sector dark matter[8, 9, 10, 11, 12]. In this scenario dark matter arises from 
a multi-component hidden sector which contains an unbroken U(l)' gauge interaction. 
The new gauge field is presumed to interact with the standard U(l)y via renormalizable 
kinetic mixing, leading to Rutherford-type elastic scattering of the dark matter particles 
off ordinary nuclei. The specific case where the hidden sector is isomorphic to the 
ordinary sector, mirror dark matter [13, 14], has also been discussed at length, and 
found to be compatible with the DAMA, CoGeNT and CRESST-II results[12] (see also 
[8, 9, 11]). More generally, though, these works indicate that there is a large class of 
such hidden sector dark matter models which are compatible with the experiments. The 
purpose of this article is to examine in detail the simplest such multi-component hidden 
sector dark matter model. It serves as a useful prototype of generic multi-component 
dark matter models of this kind. 

The outline of this article is as follows. In section 2 we review some general aspects 
of the considered hidden sector dark matter framework. In section 3 we look at the 
specific case of the model with two stable hidden sector particles, F\ and F2. In section 4 
we confront this two component hidden sector model with the latest data from DAMA, 
CoGeNT and CRESST-II experiments. In section 5, we follow up on an example near 
the combined best fit identified in section 4. In section 6, we examine the compatibility 
of this example with CDMS and XENON100 data. Finally, in section 7 we conclude. 
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2 Hidden sector dark matter 



The models we consider relegate dark matter particles to a generic hidden sector with 
an unbroken U(l)' gauge interaction. The special case where the hidden sector is 
isomorphic to the standard model sector - mirror dark matter - has been discussed in 
the context of direct detection experiments in ref.[8, 9, 11, 12]. Here we consider the 
more generic hidden sector case. This includes models with broken mirror symmetry [15] 
and many other hidden sector models[16]. More generally, hidden sector models with 
light stable particles can also be motivated from the similar inferred energy density of 
ordinary and dark matter in the Universe [1 7] . 

In this paper we consider spin-independent elastic scattering of dark matter parti- 
cles on ordinary nuclei. The dark matter particles are assumed to arise from a hidden 
sector which contains an unbroken U(l)' gauge interaction and possibly other gauge 
interactions. This means that the interactions of the theory are described by the La- 
grangian: 

£ = ^SM(e, v, u, d, B^, ...) + C dark (Fi, F 2 , .., F N , B'^, ...) + C mix (1) 

where C mix contains possible interactions connecting ordinary and hidden sector parti- 
cles. An example of Cdark, with N Dirac fermions Fj, j = 1, N, is 

Afa* = -\Kv F '^ + % ("V*" - m i) F i ( 2 ) 

where = + ig'B'^Q' is the covariant derivative and j is summed over 1, N . In 
this example the N particles Fj are each absolutely stable due to the presence of C/(l)i<g) 
^(1)2 <8> •••• <S> U (1) n global symmetries. Importantly these symmetries are not imposed; 
they are accidental, like the U(1)b and U(1)l global symmetries of the standard model. 
Evidently, having dark matter arise from a hidden sector is theoretically attractive 
because it is a natural way to introduce dark matter particles which are dark, massive, 
stable and importantly do not modify standard model physics. 

Dark matter elastic scattering depends on the cross-section, da/dEn, and on the halo 
distribution function, fp , of the dark matter species (taking the general case of N stable 
species). Both of these quantities depend crucially on the particle physics. Considering 
first the cross-section, this can have different recoil energy dependence depending on 
the type of interaction. The most common form considered in the literature is contact 
interactions which give an energy independent cross-section (excepting here the energy 
dependence due to the form factor). However, we consider Rutherford scattering which 
features non-trivial energy dependence, da/dE^ oc 1/ E\[$\. From the point of view of 
hidden sector dark matter, where the hidden sector contains an unbroken U(l)' gauge 
interaction, the latter appears to be the most natural form of interaction since it arises 
from U(1)'—U(1)y renormalizable kinetic mixing. Rutherford scattering results because 
the interaction is mediated by a massless particle (the photon). 2 

2 The case where the mediator is light is also possible, and has been studied recently in ref.[18]. 
As discussed there, in the context of the DAMA, CoGeNT and CRESST-II experiments, Rutherford 
scattering arises if the mediator has mass much less than 10 MeV, while point-like interaction arises if 
the mediator has mass much greater than 10 MeV. 



2 



The distribution function, /p., also depends on the particle physics. For dark matter 
composed of only one type of particle with mass, m, the distribution is typically assumed 
to be Maxwellian, with an effective temperature set by the galactic rotational velocity: 
T = \mv 2 ot . This leads to a distribution: / = exp(—v 2 /vl) with v = v rot . On the 
other hand, if dark matter is multi-component and self-interacting both plausible if 
dark matter arises from a hidden sector with unbroken U(l)' gauge interaction, then 
the distribution function is still Maxwellian, but the velocity dispersion depends on 
the mass of the species. Since f Fj = exp(—E/T) = exp(—^mjV 2 /T) , it follows that 

vo[Fj] = \JlT / rrij (recall the index j labels the stable particle species, j — 1, N). If 
the interactions are frequent enough so that the particles in the halo have a common 
(local) temperature, T, then it follows that the velocity dispersion depends on rrij via 
vq[Fj] oc 1/^/mJ. In fact assuming an isothermal pressure supported halo in hydrostatic 
equilibrium, one can show that [9] 



r 7-1 i iiu 

v [F j \=v rot J — (3) 




where m = J2 n j m j/ is the mean mass of the particles in the halo. Thus the 
velocity dispersion can vary greatly depending on the masses within the model. 



3 Two component hidden sector dark matter 

Consider the case where the hidden sector contains two stable U (1)' charged dark matter 
particles F\ and F 2 with masses and m F2 \F\ and F 2 can be fermionic as in the 
example Eq. (2) or alternatively bosonic] . If their abundance in the Universe arises from 
a particle-antiparticle asymmetry then two stable particles with opposite in sign, but 
not necessarily equal in magnitude, Q' charges is actually the minimal case given Q' 
neutrality of the Universe. In fact this neutrality is not just global, U{1)' electric fields 
will ensure local neutrality of a plasma containing F\ and F 2 particles. That is, 

n Fl Q' Fl + n F2 Q' F2 = (4) 

where n Fj [Q F ] is the local number density [U(l)' charge] of Fj. 

We assume that the U(l)' gauge field interacts with the standard U(l)y gauge field 
via renormalizable kinetic mixing interaction[19]: 

e ' 

r ■ = - F» v F' (5) 

mix 2cos^ ^ K) 

where F^ is the standard U(l)y gauge boson field strength tensor, and F' is the 
field strength tensor for the hidden sector U(l)'. This interaction enables the hidden 
sector U(l)' charged particles Fj to couple to ordinary photons[20] with electric charge 
9'Q' Fj e ' = e F, e - [Henceforth we define e = e F2 \. The cross-section of a F 2 particle to 
elastically scatter off an ordinary nucleus, A, presumed at rest with atomic number Z, 
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is given by [8]: 3 



da 



A 



(6) 



dE R 



EW 



where 



A = 



2ve 2 Z 2 a 2 



(7) 



m A 



and Fa^ta) is the form factor which takes into account the finite size of the nuclei. A 
simple analytic expression for the form factor, which we adopt in our numerical work, 
is the one proposed by Helm [21, 22]. 

To explain the rotation curves in spiral galaxies, we assume that the F 1 and F 2 
particles are (roughly) spherically distributed in galaxies in a pressure supported halo 
[c.f. ref.[23]]. We assume that the self-interactions are frequent enough as to thermalize 
the distributions of F\ and F 2 with a common temperature T. We further assume 
that the binding energy of atomic bound states containing Fi and F 2 particles are 
much less than this temperature, so that the F x and F 2 particles can be treated as 
two components of a plasma 4 . Typically this requires a' 2 m Fl m F2 / '{m Fl + m F2 ) keV 



As mentioned earlier, the condition of hydrostatic equilibrium relates the tempera- 
ture of the particles to the galactic rotational velocity, v rot , resulting in a mass dependent 
velocity dispersion, Eq.(3). In the two component case the mean mass of the particles 
in the galactic halo is given by 



where we have made use of Eq.(4). The interesting region of parameter space, from the 
point of view of direct detection experiments, is where m Fl <C m F2 and \Q' F \ <C \Q'f 2 \ 
so that m <C m F2 . It then follows from Eq.(3) that Vq(F 2 ) <C v^ ot . The narrow velocity 
dispersion (recall <r 2 = 3t> 2 /2) can greatly reduce the rate of F 2 interactions in higher 
threshold experiments such as XENON100[25] and CDMS/Ge [26] whilst still allowing 
F 2 to explain the signals in the lower threshold DAMA and CoGeNT experiments. 
The Fi state can be too light to be directly detected in the experiments (roughly this 
means that m Fl ~ 5 GeV), but influences the way F 2 interacts due to its effect on the 
velocity dispersion of F 2 . With these assumptions, current experiments depend on three 

3 Unless otherwise indicated, we employ natural units where fi = c = 1. 

4 The alternative case, where F\ and F2 are tightly bound into atoms has been discussed recently 
in ref. [24]. In the tightly bound limit, the interaction between ordinary matter and such atomic dark 
matter becomes point-like, and its implications for direct detection experiments are quite different 
from the case we discuss here. 



where a' = g' 2 Q' Fl Q' F2 /^n. 



n Fl m Fl + n F2 m F2 



n Fl + n F2 
Q' Fl 

m Fl - -Qf±m F2 




(8) 
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parameters, mi? 2 , m, ey £f 2 (the parameter, £p 2 w iU be defined in a moment). This two 
component hidden sector dark matter model has been briefly discussed previously in 
ref.[ll]. Our purpose here is to study it in more detail taking into account the tentative 
dark matter signal coming from the CRESST-II experiment [7] as well as the updated 
results from CoGeNT[6]. 

The differential scattering rate for Fj on a target nuclei, A, is given by 5 : 

dR AT r d ° /^( V ' V ^) | IJ3 /q\ 

N t^ Fj I -j^ '—r |v|d v (9) 



dE R 3 J\v\>v mm dE R k 

where the integration limit, v m i n , is defined by the kinematic relation: 



(m A + m Fj ) 2 E R 
\| 2m A m? F . 



(10) 



In Eq.(9), Nt is the number of target nuclei per kg of detector and n Fj is the number 
density of the halo dark matter particles Fj at the Earth's location. This number 
density can be expressed in terms of the halo mass fraction of Fj, £ Fj , and total mass 
density, p dm via n Fj = p dm £Fjm Fj (we set p dm = 0.3 GeV/cm 3 ). Also in Eq.(9) v 
denotes the velocity of the halo particles relative to the Earth and v E the velocity of 
the Earth relative to the galactic halo 6 . The halo distribution function, in the reference 
frame of the Earth, is then given by the Maxwellian distribution: 

The integral, Eq.(9), can be simplified in terms of error functions[10, 22] and solved 
numerically. 

Detector resolution effects can be incorporated by convolving the rate with a Gaus- 
sian: 

dR 1 r dR 



dE 1 ^ \phxo r 



I ^ e -{E R -E^/2al esdER (12) 

J dE R 



where E^ is the measured energy and a res defines the resolution. The unit of measured 
energy is often in keVee units (ionization/scintillation energy). For nuclear recoils, 
without channeling, keVee = keV/g, where q < 1 is the relevant quenching factor. 
Channeled events, where scattered target atoms travel down the axis and planes of the 
target crystal, have q ~ 1. In view of the theoretical study[27], we assume that the 
channelled fraction is negligible. Of course channeling may involve subtleties and the 
modeling used in ref.[27] might have missed important effects. Hence it is, of course, 
still conceivable that channeling could play an important role, which could perhaps 
modify the favored regions of parameter space. 



5 The upper limit of integration in Eq.(9) is taken as infinity since we are dealing with dark matter 
particles with potentially significant self-interactions. The self-interactions can prevent particles in the 
high velocity tail of the Maxwellian distribution from escaping the galaxy. 

6 In all numerical work we include an estimate of the Sun's peculiar velocity so that (\ve\) = v rot + 12 
km/s. 
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Bin / keV Total events Estimated background 



10.2 - 13.0 
13 - 16 
16 - 19 
19 - 25 
25 - 40 



9 

15 
11 
12 
20 



3.2 

6.1 

7.0 

11.5 

20.1 



Table 1: CRESST-II data: total number of events and estimated background. 

4 Direct detection of hidden sector dark matter 

As discussed earlier, the two component hidden sector model supposes dark matter 
arises from a hidden sector with two components F\ and F2. With mp 1 ~ 5 GeV, the 
existing direct detection experiments depend on three parameters: m F2 , e^p^, rh. It 
was shown in ref.fll] that with rh = 1 GeV and v rot = 240 km/s (fixed as a specific 
example) this model could provide a reasonable fit to the DAMA annual modulation 
signal and measured CoGeNT spectrum. Since that time, the measured CoGeNT spec- 
trum has undergone a significant 'surface event correction' [6]. Additionally, CRESST-II 
have announced results suggesting a tentative dark matter signal. Our purpose now is 
to re-examine the two component hidden sector model in light of these new experimen- 
tal developments. We shall also extend our previous analysis by considering a wider 
range of rh, v ro t values. In particular we consider rh/GeV = 1.0 and 3.0 and two values 
for the rotational velocity, v rot = 200 km/s and v rot = 240 km/s. Although v rot = 240 
km/s is representative of recent measurements of the local rotational velocity of the 
Milky Way[28] there are significant uncertainties in this quantity and it is useful to see 
how things change when v rot is varied. 

The CRESST-II experiment 

The CRESST-II experiment has collected a 730 kg-day exposure using a CaW0 4 
target [7]. The data arises from 8 detector modules, with recoil energy thresholds (keV) 
of 10.2, 12.1, 12.3, 12.9, 15.0, 15.5, 16.2, 19.0. With a view to performing a standard x 2 
analysis, we bin this data into 5 bins, summarized in table 1. This table also includes 
the expected background rate [7] (see also [29]). 

For the CRESST-II analysis, we define x 2 m the standard way: 



where Ri is the theoretically predicted rate and Bi is the estimated background in the 
jth ener gy bi n . xhe relevant rates, Ri, for the CRESST-II experiment are computed 
using Eq.(9) and Eq.(12). We = 0.3 keV[7] and assume detection efficiencies, 

e/, for the three target elements: e/ = 0.9 for O, W and e/ = 1.0 for Ca, which take 
into account their acceptance region. We have not assumed any uncertainty in the 
CRESST-II energy scale. 




(13) 
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The DAM A experiment 



The DAMA/Nal and DAMA/LIBRA experiments have accumulated data from a 
large exposure [1.17 ton-year] using a Nal target [2, 3]. Analysis of this data has yielded 
an annual modulation with phase consistent with dark matter expectations at around 
8.9o" C.L. We analyse this annual modulation signal in the energy range 2 — 8 keVee 
using 12 bins of width 0.5 keVee[3]. The detector resolution, obtained from ref.[30], has 
been included. An important issue is the nuclear recoil energy scale which is set by the 
values of the quenching factors qNa, Qi- Following our mirror dark matter analysis[12], 
we consider a range of possible values for the quenching factors: 



We minimize x 2 varying q^ ai qj over the above range of values (assumed energy inde- 
pendent for simplicity). Values of qi, q^ a outside the above range are possible and have 
been discussed in the literature. For example, the higher values q^ a ~ 0.6, qj ~ 0.3 
have been suggested by Tretyak[31]. Although the possibility of such high quenching 
factors are not specifically considered in our numerical work, qualitatively, we can say 
that higher quenching factors generally move the DAMA allowed region to lower values 
of eJ^ F2) mp 2 and also leads to interesting parameter space overlapping with that from 




CoGeNT and CRESST-II. 

The CoGeNT experiment 

The CoGeNT collaboration has recently presented data corresponding to an ex- 
posure of 0.33 x 807 kg-days in a Germanium target [6]. This update of their earlier 
exposure [5] includes an important correction for surface events which have not been 
excluded by their rise time cut [6]. We take the efficiency corrected and surface event 
corrected CoGeNT data from figure 21 of ref. [6]. We analyse this data using 15 bins of 
width 0.1 keVee in the region 0.5—2.0 keVee. We allow for a constant background, which 
we fit to the data in this energy range. We have taken into account some uncertainties 
in energy scale by minimizing the x 2 f° r CoGeNT over the variation in quenching factor, 
q Ge = 0.21 ± 0.04. 

The analysis 

We summarize the x 2 (?™n) values for the relevant data sets from each experiment 
for the two considered m values in table 2 for v rot = 200 km/s and table 3 for v rot = 240 
km/s. 

The 2a (3a) favored region of parameter space is bounded by the contours where 
X 2 {in F2 , e<Jtj? 2 ) = Xmm+4(9). In figure 1 we plot the CRESST-II 2a favored region along 
with the 3a favored region of parameter space for DAMA and CoGeNT for v ro t = 200 
km/s and with reference value of m: m — 1.0 GeV (figure la), fh = 3.0 GeV (figure 
lb). In figure 2 we repeat the exercise, but with v rot = 240 km/s. 



qNa — 



0.28 ±0.08, q! = 0.12 ±0.08. 



(14) 
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Experiment 


X 2 (min)/d.o.f. 


m/GeV 




Best fit parameters 




DAMA (annual mod.) 


5.8/10 


fh = 1.0 


rriF 2 


= 56.0 GeV, ey^F 2 = 7.2 x 10" 


-9 


CoGeNT (spectrum) 


9.4/12 


fh = 1.0 


m F2 


= 36.0 GeV, e^p~ 2 = 4.9 x 10" 


-9 


CRESST (spectrum) 


0.1/3 


fh = 1.0 


m F2 


= 80.0 GeV, e^F~ 2 = 4.9 x 10" 


-9 


DAMA (annual mod.) 


6.1/10 


fh = 3.0 


mp 2 


= 58.0 GeV, e^p~ 2 = 7.8 x 10" 


-9 


CoGeNT (spectrum) 


10.5/12 


fh = 3.0 


m F2 


= 44.0 GeV, ey^T = 5.5 x 10" 


-9 


CRESST (spectrum) 


0.3/3 


fh = 3.0 


m F2 


= 74.0 GeV, e^p~ 2 = 5.0 x 10" 


-9 



Table 2: Summary of x 2 (mm) for the relevant data sets from the DAMA, CoGeNT 
and CRESST-II experiments for two reference in values and v rot = 200 km/s. 



Experiment 


X 2 (min)/d.o.f. 


fh i 


/GeV 




Best fit parameters 




DAMA (annual mod.) 


5.0/10 


fh 


= 1.0 


m F2 


= 42.0 GeV, 




= 7.4 x 10" 


-9 


CoGeNT (spectrum) 


9.7/12 


fn 


= 1.0 


mp 2 


= 28.5 GeV, 




= 4.7 x 10" 


-9 


CRESST (spectrum) 


0.2/3 


fa 


= 1.0 


m F2 


= 57.7 GeV, 




= 4.1 x 10" 


-9 


DAMA (annual mod.) 


6.1/10 


fa 


= 3.0 


m F2 


= 46.0 GeV, 


e \JiF 2 


= 9.2 x 10" 


-9 


CoGeNT (spectrum) 


10.7/12 


fa 


= 3.0 


mp 2 


= 36.5 GeV, 




= 5.5 x 10" 


-9 


CRESST (spectrum) 


0.4/3 


fa 


= 3.0 


rriF 2 


= 39.8 GeV, 


e \[iF2 


= 4.9 x 10" 


-9 



Table 3: Summary of f (min) for the relevant data sets from the DAMA, CoGeNT 
and CRESST-II experiments for two reference fh values and v rot = 240 km/s. 
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Figure la: DAMA, CoGeNT and CRESST-II favored regions of parameter space in the two 
component hidden sector model. The reference point v ro t = 200 km/s and m = 1.0 GeV is 
assumed. 
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Figure lb: Same as figure la, except that m = 3.0 GeV is assumed. 
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Figure 2a: DAMA, CoGeNT and CRESST-II favored regions of parameter space in the two 
component hidden sector model. The reference point v ro t = 240 km/s and m = 1.0 GeV is 
assumed. 
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Figure 2b: Same as figure 2a, except that m = 3.0 GeV is assumed. 



10 



The figures demonstrate that there are quite large regions of parameter space where 
each experiment can be explained within this hidden sector framework. For DAMA, 
the signal is dominated by F 2 — Na scattering if ~ 40 GeV, while if mp 2 ~ 40 
GeV then both F 2 — Na and F 2 — I scattering conspire to produce the signal. Also, the 
signal arises from scattering of target nuclei with F 2 particles from the body of their 
Maxwellian halo distribution (rather than, say, the tail) and thus is reasonably stable, 
and shows little dependence on the velocity dispersion i.e. on m. 

For CoGeNT, the spectrum is consistent with dR/dE R oc l/E^, predicted from the 
energy dependence of the Rutherford cross section. This is the reason why the CoGeNT 
spectrum is reproduced for a large range of F 2 mass. In other words, the shape of the 
CoGeNT spectrum arises from dynamics rather than kinematic effects in this model. 
Again there is little dependence on the m parameter. 

Finally, in the case of CRESST-II, the dominant signal contribution arises from 
F 2 — Ca scattering for m F2 ~ 50 GeV. For such masses, F 2 particles in the body of their 
Maxwellian halo distribution can scatter to produce nuclear recoils above the 10 keV 
threshold. As the mass of F 2 is lowered, mp 2 ~ 40 — 50 GeV, a recoil above threshold 
can only occur for F 2 particles in the tail of the Maxwellian halo distribution. For 
this reason, there is quite a bit of dependence of the CRESST-II allowed region on the 
velocity dispersion, i.e. rh in this mass range. 

The figures show that the allowed region of parameters favoured by DAMA, CoGeNT 
and CRESST-II are very similar with some degree of overlap. This occurs ignoring a 
variety of possible systematic effects, such as the possibility of channeled events in 
DAMA (which can lower the DAMA favoured values of e^£p 2 ), and uncertainties in 
the surface event correction factor for CoGeNT (which could raise or lower the CoGeNT 
favoured values of ey^^"). Thus realistically there is quite a large region of parameter 
space that is possible, much larger than the overlapping allowed region, even assuming 
all three experiments have detected dark matter interactions. It is difficult, though, to 
quantify all the possible systematic effects and do an exhaustive analysis. Instead, we 
shall hope to gain some insight by examining a particular point in parameter space. 

5 An example near the combined best fit of DAMA, 
CoGeNT and CRESST-II 

In view of potential systematic uncertainties between the experiments we shall resist 
the temptation to fit the combined DAMA, CoGeNT and CRESST-II data. Instead 
we shall consider an example reference point PI located near the overlapping allowed 
regions indicated in figure la. 

PI : m F2 = 50 GeV, = 5.7 x 10~ 9 , rh = 1.0 GeV, v rot = 200 km/s (15) 

The DAMA annual modulation signal for this example is given in figure 3a, the CoGeNT 
spectrum in figure 3b, and CRESST-II spectrum in figure 3c. 
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Figure 3a: DAMA annual modulation spectrum for two component hidden sector dark matter 
with parameters PI (solid line). The separate contributions from scattering off sodium (dotted 
line) and Iodine (dashed-dotted line) are also shown. In this example q^ a = 0.36, qj = 0.20. 
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Figure 3b: CoGeNT spectrum for two component hidden sector dark matter with parameters 
PI. In this example qc e = 0.17. 
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Figure 3c: CRESST spectrum for two component hidden sector dark matter for the example 
PI. 



Figures 3 clearly demonstrates that two component hidden sector dark matter can 
simultaneously explain DAMA, CoGeNT and CRESST-II data. The apparent bump 
in the DAMA data can be reproduced extremely well and the minimal model predicts 
a change in sign of the DAMA annual modulation signal at low energies. Note though, 
that if there are additional particles, F$, of intermediate mass ~ 20 GeV, then their 
positive contribution to the annual modulation can outweigh the negative contribution 
from F 2 . An example, given in the context of mirror dark matter with F 2 ~ Fe' and 
F 3 ~ O', was discussed in ref.[12]. Figure 3b indicates that the shape of the CoGeNT 
spectrum is nicely consistent with the dR/dE R oc l/E^ predicted from the energy 
dependence of the Rutherford cross section. The event rate is predicted to continue to 
rise as ~ 1/-E# as the recoil energy is reduced below the current CoGeNT threshold, 
until the threshold of a lighter component is reached, whereby the rate can jump even 
higher. These effects can be probed by TEXONO, C-4 and possibly other experiments. 

Future data from DAMA, CoGeNT, CRESST-II and other experiments will obvi- 
ously be able to constrain the parameter space within this hidden sector framework. 
As discussed recently [32], a particularly striking diurnal modulation signal should be 
observable for a detector located in the southern hemisphere, and perhaps even in 
the northern hemisphere at low latitudes, such as detectors in Jin-Ping underground 
laboratory. In the meantime, we must rely on annual modulation and spectrum data. 

In figure 4 we give the total predicted spectrum for DAMA/LIBRA for the reference 
point, PI. Figure 4 also shows the 'single hit' event rate recorded in DAMA/LIBRA [3]. 
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The figure indicates that the sharp rise in the predicted dark matter interaction rate 
could potentially be differentiated from the background if the DAMA threshold is low- 
ered below 2 keVee. In figure 5 we show the predicted annual modulation spectrum 
for CoGeNT for the same reference point, PI. Clearly, the initial annual modulation 
amplitude measured by CoGeNT to be A « 0.46 ±0.17 cpd/kg/keVee averaged over the 
energy range: 0.5 < -E(keVee) < 3.0, is much larger than that predicted by our example 
point. Also, we find that the annual modulation changes sign at low energies. This 
feature is not supported by CoGeNT's initial measurement. However the energy where 
the modulation changes sign can be reduced if the mass of F 2 is lowered, which can 
alleviate this discrepancy. Alternatively, there can be an additional dark component, 
i*3, of mass ~ 20 GeV. As discussed earlier in the context of DAMA, such a component 
can give a positive contribution to the annual modulation which can outweigh the neg- 
ative contribution from F 2 . Clearly future measurements of the annual modulation by 
CoGeNT, C-4, CDEX, and other experiments will be of crucial importance. 
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Figure 4: DAMA spectrum for hidden sector dark matter with parameters PI. In this example 
q Na = 0.36, QI = 0.20. 
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Figure 5: CoGeNT annual modulation spectrum for hidden sector dark matter with parame- 
ters PI. In this example qc e = 0.17. 

6 Constraints from CDMS and XENON100 

The CDMS/Ge[26], CDMS/Si[33] and XENON100[25] experiments have reported null 
results in their dark matter search. Here, we examine the compatibility of the considered 
hidden sector dark matter model with these null results 7 . The constraints from the 
CDMS and XENON100 experiments depend sensitively on the recoil energy threshold 
which (typically) has at least 20% uncertainty. The XENON100 experiment, which 
has the largest exposure, has particularly poor knowledge of energy calibration in the 
important low energy region (see ref.[38, 39] for some relevant discussions). Given this 
situation, we examine the compatibility issue by estimating the energy threshold for 
which the parameter point PI can be excluded at 95% C.L. We have taken into account 
the relevant detection efficiencies, exposure time etc, although we conservatively ignore 
effects of detector resolution for the XENON 100 experiment. Our results are given in 
table 4. Note that the 95% C.L. limits in the second column arise since CDMS/Ge, 
CDMS/Si and XENON100 observed 2, 0, 2 events respectively in the low energy region. 



7 There are also lower threshold analysis by the XENON10[34] and CDMS collaborations [35]. How- 
ever it has been argued [36] that neither analysis can exclude light dark matter (and by extension, 
hidden sector dark matter considered here, which has similar event rates at low energies) when sys- 
tematic uncertainties are properly taken into account. Interestingly, a recent analysis[37] has found 
that the low energy CDMS data are actually fully consistent with CoGeNT's observed low energy 
excess rate, adding weight to the dark matter interpretation of this excess. 
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Experiment 


95% C.L limit 


Eg threshold 


E nom threshold 


100 


E nom 
E R 




CDMS/Ge 


6.3 


12.8 keV 


10.0 keV 




22 % 


CDMS/Si 


3.0 


8.8 keV 


7.0 keV 




20 % 


XENON100 


6.3 


13.0 keV 


6.6 keV 




49 % 



Table 4: Estimated energy threshold (Eg) for which the CDMS/Ge, CDMS/Si and 
XENON100 experiments are consistent with the hidden sector dark matter expecta- 
tions assuming the example point, PI. E™ m is the nominal energy threshold of each 
experiment. The last column gives the percentage at which the energy threshold needs 
to have been underestimated for the point PI to be consistent with these experiments 
at 95% C.L. 

Table 4 indicates that hidden sector dark matter with parameters given by the ex- 
ample PI is consistent with the CDMS/Ge, CDMS/Si [XENON100] data provided that 
the energy scale has been underestimated by of order 20% [50%]. This indicates that 
the hidden sector dark matter explanation of the DAMA, CoGeNT and CRESST-II 
experiments has some level of tension - especially with the XENON100 results. Our ex- 
ample point assumed fh = 1.0 GeV. Raising fh will increase the tension of hidden sector 
dark matter with the null results of CDMS and XENON100 experiments. Lowering fh 
can improve the situation, but only moderately so, since fh = 1 GeV is already very 
low. Although PI was just an example point, it seems that some level of tension exists 
between the null results of XENON100 and the hidden sector dark matter explanation 
of the DAMA, CoGeNT and CRESST-II experiments. 

One can envisage several possible ways in which this tension might be alleviated. 
For example, it is possible that there is an issue with the calibration of the XENON100 
apparatus. Determining the recoil energy scale in the XENON100 detector is nontrivial 
and it seems possible that this scale might be have a factor ~ 2 uncertainty [38, 39]. 
Another possibility is that F 2 has a somewhat lower mass than given in our example. 
This option would be especially relevant if we were to ignore the CRESST-II excess. 
For example the parameter point: 

P2 : m F2 = 40 GeV, ey^T = 5.7 x 10~ 9 , fh = 1.0 GeV, v rot = 190 km/s (16) 

yields a ^(dama) = 18.9 for 12 data points and ^(cogent) = 22.8 for 15 data points, 
but gives only around 1 event for the CRESST-II exposure. This is a reasonable fit for 
DAMA and CoGeNT, considering that the shape of the distributions are fit well, only 
the overall normalization is not [CoGeNT (DAMA) prefers slightly smaller (larger) 

e \J^F 2 ) with e^J~£jr 2 = 5.7 x 10~ 9 a compromise]. The different normalizations might 
easily be due to systematic effects not included in our analysis. In table 5 we estimate 
the energy threshold for which the CDMS/Ge, CDMS/Si and XENON100 experiments 
are consistent with hidden sector dark matter expectations assuming the example point 
PI. 

Finally, the KIMS experiment [40] is currently operating with a ~ 100 kg Csl target. 
Using a Pulse Shape Analysis (PSA), they have announced a (2a) limit on the rate of 



16 



Experiment 


95% C.L limit 


Eg threshold 


E nom threshold 


100 


E nom 
E R 




CDMS/Ge 


6.3 


9.8 keV 


10.0 keV 




-2% 


CDMS/Si 


3.0 


7.3 keV 


7.0 keV 




4% 


XENON100 


6.3 


9.4 keV 


6.6 keV 




30% 



Table 5: Same as Table 4, except for the parameter point P2. 



dark matter interactions: R < 0.02 cpd/kg/keVee for 3 < ^(keVee) < 4. We find that 
the predicted rate of Iodine interactions in a Csl target, for the parameter point PI, 
is R ~ 0.007 cpd/kg/keVee in the 3 < En(keVee) < 4 bin and falls sharply at higher 
energies. [Assuming here the resolution and qj quenching factor for KIMS are the same 
as in the DAMA/LIBRA experiment]. Thus, their PSA does not currently constrain 
hidden sector dark matter. However, the KIMS experiment also hopes to measure or 
constrain the annual modulation amplitude. These results are expected to be reported 
shortly. We note here that KIMS can potentially see a positive signal for a significant 
portion of parameter space (c.f. dashed-dotted line in figure 3a). 

7 Conclusion 

We have examined the data from the DAMA, CoGeNT and CRESST-II experiments 
in the context of multi-component hidden sector dark matter models. The models 
considered feature a hidden sector with two or more stable particles charged under an 
unbroken U(l)' gauge interaction. The new gauge field can interact with the standard 
U(l)y via renormalizable kinetic mixing, leading to Rutherford-type elastic scattering 
of the dark matter particles off ordinary nuclei. We considered the simplest generic 
model of this type, with a hidden sector composed of two stable particles, F 1 and F 2 . 

We have found that the two component hidden sector dark matter model can si- 
multaneously explain the DAMA, CoGeNT and CRESST-II data. This explanation 
is consistent with the null results of the CDMS/Si and CDMS/Ge experiments when 
~ 20% systematic uncertainties in energy scale are considered. There is some tension 
with the XENON100 experiment. The favoured parameter regions are typically consis- 
tent with the most recent XENON100 results only if the XENON100 energy threshold 
is around a factor of two higher than given by the XENON100 collaboration. 
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